Variational solution of the Yang-Mills Schrodinger equation in Coulomb gauge* 
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The Yang-Mills Schrodinger equation is solved in Coulomb gauge for the vacuum by the varia- 
tional principle using an ansatz for the wave functional, which is strongly peaked at the Gribov 
horizon. A coupled set of Schwinger-Dyson equations for the gluon and ghost propagators in the 
Yang-Mills vacuum as well as for the curvature of gauge orbit space is derived and solved in one- 
loop approximation. We find an infrared suppressed gluon propagator, an infrared singular ghost 
propagator and a almost linearly rising confinement potential. 
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I. INTRODUCTION 



To understand the low energy sector of QCD is one 
of the most challenging problems in quantum field 
theory. Nowadays, quantum field theory and, in par- 
ticular, QCD is usually studied within the functional 
integral approach. This approach is advantagous for a 
perturbative calculation, where it leads automatically 
to a Feynman diagrammatic expansion. Within this 
approach asymptotic freedom of QCD was shown p|, 
which manisfests itself in deep inelastic scattering ex- 
periments. In addition, the functional integral approach 
is the basis for numerical lattice calculations 0, the 
only rigorous non-perturbative approach available at 
the moment. These lattice methods have provided 
considerable insights into the nature of the Yang-Mills 
vacuum. Lattice investigations performed over the last 
decade, have accumulated evidence for two confinement 
scenarios: the dual Meissner effect based on a condensate 
of magnetic monopoles and the vortex condensation 
picture 4] (For a recent review see ref. Q)- In both 
cases the Yang-Mills functional integral is dominated 
in the infrared sector by topological field configurations 
(magnetic monopoles |(j or center vortices which 
seem to account for the string tension, i.e. for the 
confining force. Yet another confinement mechanism 
was proposed by Gribov further elaborated by 
Zwanziger || and tested in lattice calculations 0| . This 
mechanism is based on the infrared dominance of the 
field configurations near the Gribov horizon in Coulomb 
gauge. This mechanism of confinement is compatible 
with the center vortex and magnetic pictures, given the 
fact, that lattice center vortex and magnetic monopolc 
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configurations lie on the Gribov horizon 
Despite the great successes of lattice calculations in the 
exploration of strong interaction physics [T^ | , a complete 
understanding of the Yang-Mills theory will probably 
not be provided by the lattice simulations alone, but 
requires also analytical tools. Despite of its success 
in quantum field theory, in particular in perturbation 
theory and its appeal in semi-classical and topological 
considerations |l3j |. the path integral approach may 
not be the most economic method for analytic studies 
of non-perturbative physics. As an example consider 
the hydrogen atom. Calculating its electron spectrum 
exactly in the path integral approach is exceedingly com- 
plicated [Til ], while the exact spectrum can be obtained 
easily by solving the Schrodinger equation. One might 
therefore wonder, whether the Schrodinger equation is 
also the appropriate tool to study the low-energy sector 
of Yang-Mills theory, and in particular of QCD. 
The Yang-Mills Schrodinger equation is based on the 
canonical quantization in Weyl gauge Aq = [l5| . 
where Gaufi' law has to be enforced as a constraint 
to the wave functional to guarantee gauge invariance. 
The implementation of Gaufi' law is crucial. This is 
because any violation of Gaufi' law generates spurious 
color charges during the time evolution. These spurious 
color charges can screen the actual color charges and 
thereby spoil confinement 36]. Several approaches 
have been advocated to explicitly resolve Gaufi' law by 
changing variables resulting in a description in terms 
of a reduced number of unconstrained variables. This 
can be accomplished either by choosing a priori gauge 
invariant variables |l7j or by fixing the gauge, for 
example, to unitary gauge 0] , to Coulomb gauge 0] or 
to a modified version of axial gauge [2j}. In particular, 
the Yang-Mills Hamiltonian resulting after eliminating 
the gauge degrees of freedom in Coulomb gauge, was 
derived in ref. [l^. Alternatively, one has attempted 
to project the Yang-Mills wave functional onto gauge 
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invariant states (which a priori fulfill Gaufi' law) |21j . 
The equivalence between gauge fixing and projection 
onto gauge invariant states can be seen by noticing 
that the Faddeev-Popov determinant provides the Haar 
measure of the gauge group 22]. 

In this paper we will variationally solve the stationary 
Yang-Mills Schrodinger equation in Coulomb gauge for 
the vacuum. Such a variational approach was previously 
studied in refs. [2j|, |24|, where a Gaussian ansatz for 
the Yang-Mills wave functional was used. Our approach 
is co ncep tually simuiar to, but differs essentially from 
refs. [23,|2J] in two respects: i) we use a different ansatz 
for the trial wave functional and ii) we include fully 
the curvature of the space of gauge orbits induced by 
the Faddeev-Popov determinant. We use a vacuum 
wave functional, which is strongly peaked at the Gribov 
horizon. Such a wave functional is motivated by the 
results of ref. [26j and by the fact, that the dominant 
infrared degrees of freedom like center vortices lie on the 
Gribov horizon ^lj. The Faddeev-Popov determinant 
was completely i gno red in ref. |23| and only partly 
included in ref. |24j. (In addition, ref. [24| uses the 
angular approximation). We will find however, that a 
full inclusion of the curvature induced by the Faddeev- 
Popov determinant is absolutely crucial and vital for the 
infrared regime and, in particular, for the confinement 
property of the Yang-Mills theory. The organization of 
the paper is as follows: 

In section [H] we briefly review the hamiltonian formula- 
tion of Yang-Mills theory in Coulomb gauge and fix our 
notation. In section lTTTl we present our vacuum wave func- 
tional and calculate the relevant expectation values. The 
vacuum energy functional is calculated and minimized in 
section |W] resulting in a set of coupled Schwinger-Dyson 
equations for the gluon energy, the ghost and Coulomb 
form factors and the curvature in gauge orbit space. The 
asymptotic behaviours of these quantities in both the 
ultraviolet and infrared regimes are investigated in sec- 
tion Section IVII is devoted to the renormalization of 
the Schwinger-Dyson equations. The numerical solutions 
to the renormalizcd Schwinger-Dyson equations are pre- 
sented in section IVIII Finally in section IVIIII we present 
our results for the static Coulomb potential. Our conclu- 
sions are given in section liXJ A short summary of our 
results has been previously reported in ref. |16| . 



II. HAMILTONIAN FORMULATION OF 
YANG-MILLS THEORY IN COULOMB GAUGE 



The canonical quantization of gauge theory is usually 
performed in the Weyl gauge A = 0. In this gauge 
the spatial components of the gauge field A(x) are the 
"cartesian" coordinates and the corresponding canoni- 



cally conjugated momenta 

defined by the equal time commutation relation 

[^(x),n?(y)] =i6 ab 5 kl S(x-y) (2) 

represent the color electric field. The Yang-Mills Hamil- 
tonian is then given by 

H = \J ^[n«(x) 2 +S«(x) 2 ] , (3) 

where 

B k = ^' ,,,,/•,,. , gFij = [Di,Dj] , D t =d l + gA l (4) 

is the colormagnetic field with Fij being the non-Abelian 
field strength and Dj the covariant derivative. We use 
anti-hermitian generators T a of the gauge group with 
normalization tr(T a T b ) = ~\5 ab {A = A a T a ). 
The Hamiltonian J3J is invariant under spatial gauge 
transformations Z7(x) 

A -» A u = UAU^ + -UdU ] . (5) 
9 

Accordingly the Yang-Mills wave functional ^f[A] = 
(A\^f) can change only by a phase, which is given by 

y[A u ] =e ien ^*[A] , (6) 

where O is the vacuum angle and n [U] denotes the wind- 
ing number of the mapping E/(x) from the 3-dimensional 
space R 3 (compactified to S3) into the gauge group. In 
this paper we will not be concerned with topological as- 
pects of gauge theory and confine ourselves to "small" 
gauge transformations with n[U] = 0. 
Invariancc of the wave functional under "small" gauge 
transformations is guaranteed by Gaufi' law 

D k (A)IL k V[A] = Pm , (7) 

where p m is the color density of the matter (quark) field 
and Di(A) denotes the covariant derivative in the adjoint 
representation of the gauge group 

Di(A) = di + gAi , Ai - A$f a , (f a ) bc = f bac (8) 

with f abc being the structure constant of the gauge 
group. Throughout the paper we use a hat """ to de- 
note quantities defined in the adjoint representation of 
the gauge group. The Gaufi' law constraint Q on the 
Yang-Mills wave functional can be resolved by fixing the 
residual gauge invariancc, cq. Q. This eliminates the 
unphysical gauge degrees of freedom and amounts to a 
change of coordinates from the "cartesian" coordinates 
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to curvilinear ones, which introduces a non-trivial Jaco- 
bian (Faddeev-Popov determinant) »7[A]. In this paper 
we shall use the Coulomb gauge 

d k A k (x) = . (9) 

In this gauge the physical degrees of freedom are the 
transversal components of the gauge fields 

A x (x)=t ij ( X )A j (x), (10) 

where 

Mx) = %-^ (ii) 



H 



where 

J[A X ] = Det (-di A) (15) 

is the Faddeev-Popov determinant with Dk = Dk(A x ) . 
The kinetic part of the hamiltonian (first term) is remi- 
niscent to the functional version of the Laplace-Beltrami 
operator in curvilinear space. The magnetic energy (sec- 
ond term) represents the potential for the gauge field. Fi- 
nally, the last term arises by expressing the kinetic part 
of eq. ® in terms of the transversal gauge potentials 
(curvilinear coordinates) satisfying the Coulomb gauge, 
i.e. by resolving GauB' law |l9|. This term describes the 
interaction of static non-Abelian color (electric) charges 
with density 

p°(x)=-i i L " i (x)Il*(x)+p^(x) (16) 

through the Coulomb propagator 

F ab (x,x') = (xa| (-D^y 1 (~3 2 ) (-D&y 1 |x'6)(17) 

which is the non-Abelian counter part of the usual 
Coulomb propagator (x|^ip-|x') in QED and reduces to 

the latter for the perturbative vacuum A 1 - = 0. If not 
stated otherwise, in the following we will assume p m = 0. 

Since the Hamiltonian l|14l) does not change when J^A^] 
is multiplied by a constant we can rescale the Jacobian 
J by the irrelevant constant Det (— d 2 ) 

Det (-Didi) 
^ = De\ { -V) ^ 



is the transversal projector. The corresponding canoni- 
cal momentum is given by 

*"w-«*w«4r3afe (12) 

and satisfies the equal-time commutation relation 

[A Xa (x) , n| b (x')] = i5 a %(x)S (x - x') (13) 

In Coulomb gauge © the Yang-Mills Hamiltonian is 
given by 0] 

(14) 



such that J[A X = 0] = 1. 

The Faddeev-Popov matrix ( — Did^ represents the met- 
ric tensor in the color space of gauge connections sat- 
isfying the Coulomb gauge. Accordingly its determi- 
nant enters the measure of the integration over the space 
of transversal gauge connections and the matrix ele- 
ment of an observable 0[A X , H] between wave functionals 
fi^], V 2 [A X ] is defined by 

{*i\0\9 2 ) = J VA x J[A x M[A x ]0[A x ,H]^ 2 [A x ] . 

(19) 

Here, the integration is over transversal gauge potentials 
A 1 - only. The same expression (|19H is obtained by 
starting from the matrix element in Weyl gauge (with 
gauge invariant wave functionals) and implementing the 
Coulomb gauge by means of the Faddeev-Popov method. 

Like in the treatment of spherically symmetric systems in 
quantum mechanics it is convenient to introduce "radial" 
wave functions by 

$[A X ] = J^[A X ]$[A X ] (20) 

and transform accordingly the observables 

d = ji[A x ]OJ-i[A :i } . (21) 

Formally, the Jacobian then disappears from the matrix 
element 

(*i|0|* 2 ) = / DA X *1[A X ] 6 V 2 [A X ] := (*i|0|* 2 ) -(22) 



~ J d 3 x [j-\A x ]U«{*)J[A x ]Ti<t{K)+B«{Kf} 
+ y J d 3 x J (fx'J-^A^p^F^ (x,x') J[A x ]p b (x') 
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In the present paper we are interested in the vacuum 
structure of Yang-Mills theory. The vacuum wave func- 
tional, defined as solution to the Yang-Mills Schrodinger 
equation 



(23) 



To work out the expectation value of the Hamiltonian 
H = J^HJ~? it is convenient to perform a par- 
tial (functional) integration in the kinetic and Coulomb 
terms. Assuming as usual that the emerging surface 
terms vanish, one finds 



for the lowest energy eigenstate, can be obtained from 
the variational principle 



(tf |JT|¥) = = E k + E P + E C 



(25) 



(*!*) 



(24) where 



E k = - I PA 



d 3 3 



ut(-x)^[A- L ] n?(x)#L4 



E p = ^ I VA 1 ^ I d 3 x ^[A^Bfixf^lA^ 



VA J 



d 3 x / d 3 x' 



n?(x)*[A x ] A ±ca (x)F oh (x,x')A ±M (x') n^x')*^- 1 ] 



(26) 
(27) 
(28) 



where 14 is defined below in eq. (|29JI . Note, integration 
is over the transversal part A 1 - only. In order to pre- 
vent the equations from getting cluttered we will in the 
following often write A instead of A 1 -, and 14 instead of 
IT 1 , but it will be always clear from the context, when 
the transversal part is meant. 

The transformed momentum operator 14° (x) is explicitly 
given by 

fL a k (x)=J l *[A x }ILi a (x.)J-3[A x ] 

= n^(x)-in ±a (x)inj[^] 

= I4^(x) + |i H (x) tr [f a (3fG(y,x)) y=x , 

(29) 

where we have introduced the inverse of the Faddeev- 
Popov operator 

G=( K -D l d l )~ 1 , (30) 
which is a matrix in color and coordinate space 

(xa|G|x'6) = (xa| (-Aft) 1 |x'6) := G ab (x,x') . 

Its vacuum expectation value (to be defined below) rep- 
resents the ghost propagator. Expanding this quantity 
in terms of the gauge field, we obtain 

G = (-d 2 - gAidi) =GoY, (dAidiGo) , (31) 

n=0 

where 

(32) 



is the free ghost propagator, which is nothing but the 
ordinary static Coulomb propagator in QED. 
From the expansion, eq. I|31|) . it is seen, that the ghost 
propagator satisfies the identity 



G = Go + GogAAG , 



(33) 



where we have used the usual short hand matrix notation 
in functional and color space. 



III. THE VACUUM WAVE FUNCTIONAL AND 
PROPAGATORS 



So far all our considerations have been exact, i.e. no 
approximation has been introduced. To proceed further, 
we have to specify the (vacuum) wave functional SSf[A], 



G (x,x') = (x|(-9 2 ) 1 



Am x — x' 



The vacuum wave functional 



Inspired by the known exact wave functional of QED in 
Coulomb gauge |25|. we will consider a Gaussian ansatz 
in the transversal gauge fields 
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V[A] = (A\uj) = Wexp 



d 3 x'Ai a (x)cj{x,x')A± a (x') 



(34) 



where J\f is a normalization constant chosen so that 
(^1^) = 1. By translational and rotational invariancc 
the integral kernel in the exponent of the Gaussian wave 
functional ui(x, x') can depend only on |x — x'|. For sim- 
plicity we have chosen the integral kernel oj(x, x') to be 
a color and Lorentz scalar. This is justified by isotropy 
of color and Lorentz space. 

Of course, such a Gaussian ansatz is ad hoc at this stage 
and can be justified only a posteriori. Let us also stress, 
that we are using the Gaussian ansatz, eq. (|34f> . for the 
radial wave function ^[A], eq. (|20() . which is normal- 
ized with a "flat" integration measure, see eq. 1|22|1 . The 
original wave function 

*[A] = J~^[A]9[A] (35) 

contains besides the Gaussian *ff[A] (|34() an infinite 
power series in the gauge potential A. Furthermore, 
since the Jacobian J[A] (Faddeev-Popov determinant) 
vanishes at the Gribov horizon, our wave functional is 
strongly peaked at the Gribov horizon. It is well known, 
that the infrared dominant configurations come precisely 
from the Gribov horizon [2y|. In addition, the center 
vortices, which are believed to be the "confiner" in the 
Yang-Mills vacuum Q all live on the Gribov horizon . 
Furthermore, the wave functional (|35|l . being divergent 
on the Gribov horizon, identifies all gauge configurations 
on the Gribov horizon, in particular those which are 
gauge copies of the same orbit. This identification is 
absolutely necessary to preserve gauge invariance. (The 
idcntificaton of all gauge configurations on the Gribov 
horizon also topologically compactifies the first Gribov 
region.) Therefore we prefer to make the Gaussian 
ansatz for the radial wave function *&[A], instead for 
*&[A}. Our wave functional thus drastically differs from 
the one used in refs. |2^| , [24| , where an Gaussian ansatz 



was used for ^f[A]. 

We should also mention that the wave functional al- 
though being divergent at the Gribov horizon it is obvi- 
ously normalizable. In principle, a wave functional being 
peaked at the Gribov horizon can of course have a more 
general form than the one given by eq. (|35fl . A somewhat 
more general wave functional 

=J- a [A]*[A] (36) 

would leave the power a as a variational parameter, 
which is then determined by minimizing the vacuum en- 
ergy (density). This would lead to a more optimized wave 
functional, which for a > still expresses the dominance 
of the field configurations on the Gribov horizon. Such 
investigations are under way. In the present paper we 
restrict however ourselves to a = i, which simplifies the 
calculations a lot. 

We will use the Gaussian ansatz l|34(l as a trial wave func- 
tion for the Yang-Mills vacuum and determine the inte- 
gral kernel w(x, x') from the variational principle, mini- 
mizing the vacuum energy density. The use of the Gaus- 
sian wave functional makes the calculation feasible in the 
sense that Wick's theorem holds: the expectation value 
of an ensemble of field operators can be expressed by the 
free (static) gluon propagator 

(4^(x)A+V)>. := H4^(xM| b (x')M 

= ^%(x) W - 1 (x,x'). (37) 

The expectation value of an odd number of field operators 
obviously vanishes for the Gaussian wave functional. 
To facilitate the evaluation of expectation values of field 
operators, we introduce the generating functional 



Z\j] = (§|exp 



I*) 



W 2 J DA ± exp - J (fx J dV^ a (x)w(x,x')A^°(x') + J < 



g^jT(x)^(x) 



(38) 



where the normalization constant TV guarantees that obtains 
Z\j = 0] = 1- Carrying out the Gaussian integral one 
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Z[j] = exp 



d 3 X / dV#(x)ty(x)w- 1 (x,X / )#( X ') 



(39) 



r 



The expectation value of any functional of the gauge field, in the following, can be obtained by using the identity 
0[A], is then given by 



(0[A]) U = (u\0[A]\u) = (0 



m ) ■ (40) 



(41) 



y=o 



This is basically the functional form of Wick's theorem. which can be proved by Fourier transformation. Apply- 
An alternative form of this theorem, which will be uscfull ing this to eq. (|40|) we obtain 



(0[A]) U 



exp 



d 3 x / d 3 x' 



^w^ (x) ^ I(x ' x,) 



SAf a (x') 



0[A] 



(42) 



A=0 



r 



Since the kernel w(x, x') depends only on |x — x'|, it is 
convenient to go to momentum space by Fourier trans- 
formation 



,4 a (x) = 



d 3 k 
(2^)3 



vkx 



A a (k) 



(43) 



In momentum space the wave functional l|34(l reads 
(A\lj) = A/" exp ' 1 ' ] 1 



, ^(k)^-(k)a;(k)^(-k) 



(44) 



and the free gluon propagator l|37() becomes 



(u;\At(k)A J (-k')\u J )=6 



a6 ^(2 7 r) 3 5(k-k') ) (45) 



where w(k) represents the energy of a gluon with momen- 
tum k. In order that our trial wave function, eq. (|34ll . is 
integrable w(k) has to be strictly positive definite. 



B. The ghost propagator and the ghost-gluon 
vertex 



which we refer to as ghost propagator, although we will 
not explicitly introduce ghost fields. They are not needed 
in the present operator approach. In the functional in- 
tegral in Coulomb gauge, eq. would enter as the 
propagator of the ghost field. Taking the expectation 
value of eq. (|33|l and using the fact, that the free ghost 
propagator Go does not depend on the dynamical gauge 
field, so that Go\u>) = \oj)Gq, we obtain 



G u = Go + G (gA i d i G 



(47) 



The latter expectation value can be worked out, in prin- 
ciple, by inserting the expansion (|31|l and using Wick's 
theorem. The result can be put into a compact form by 
defining the one particle irreducible ghost self-energy E 
by 



[gAidiGj^ = E(G) W = EG W . (48) 
It is then seen that G u satisfies the usual Dyson equation 



For later use, let us consider the vacuum expectation 
value of the inverse Faddeev-Popov operator 130|) 



(46) 



Gu; — Gn + GnSG^ 



(49) 



We will later also need the functional derivative of the 
ghost propagator 



M£(x) 



G(xi,x 2 ; 



dV /d 3 y2 G(x 1 ,y 1 )^|^G(y 2 ,x 2 ) 



(50) 
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FIG. 1: Diagrammatic expansion of the ghost-gluon vertex. 
Throughout the paper full and curly lines stand, respectively, 
for the full ghost and gluon propagators. Furthermore, dots 
and fat dots represent, repectively, bare and full ghost-gluon 
vertices. 



This expectation value can in principle be evaluated by 
using the expansion Q49JI for the ghost propagator and 
applying Wick's theorem. The emerging Feynman dia- 
grams have the generic structure illustrated in figure ^ 
They describe the propagation of the ghost followed by an 
interaction with an external gluon and subsequent prop- 
agation of the ghost. We therefore define the ghost-gluon 
vertex T£(x) by 



M°(x) G ^ 1|]C ^) = / d3yi J d3 ^G , a ) (xi,yi)r^(y 1 ,y 2 ;x)G w (y 2 ,x 2 ) 



(51) 



The ghost-gluon vertex T^(x) is given by the series of which we refer to as bare vertex r°' a (x), is given by 
diagrams shown in figure ^ Comparison of eqs. (|50fl and 
(|51|l shows, that the leading order contribution to (x) , 



= g^(y)'5(y-x 1 )f a 5f 1 ( 5(x 1 -x 2 ) (52) 

or in momentum space 

iTCxr^y) = / t^/ |^ r r^ k K k(y - Xl) • eiq(Xl " X2) (53) 

by 

r°' Q (q,k) = igf a t kl (k) qi . (54) 



C. The ghost self-energy 



The ghost self-energy £ is conveniently evaluated from its defining equation l|48|l by using Wick's theorem in the form 
ofeq. 



SG W = (gA k d k G) = {exp 



J/ A / A '^ (X)W " (X ' X ' } 6A^j 



gAAG 



| exp (ijd^Jd^-^^Uj^ (x,x') jjn^jj >° AA 



A=0 

a] 



(55) 



A=0 



Using now the relation 



and expressing j^gA k d k via eq. (|3U|) as — j^G 1 , where 



=/(A) 



(56) 



the latter quantity represents in view of eq. the bare 



ghost-gluon vertex r° , we find 



J d 3 x 3 j:(x 1 ,x 3 )G LU {x 3 ,x 2 ) = d 3 y 1 y"d 3 ?/ 2 ^ rf 3 Z3 M 1 ^^) ^'^ 1 ^ * ( yi ' y2 ) r °' C ( Xl>X3 ' y2 ) 



exp 



\ [ d 3 x f d 3 x' — ^— - tiAx)u~ l (x.x'l — r^- 



G(x 3 ,x 2 U .(57) 

J A=0 



Since T° and a; are independent of A the variational derivatives act only on G(x 3 ,x 2 ) and we obtain 

J d 3 x 3 S(xi,X3)G w (x 3 ,x 2 ) = \ J d 3 yi J d 3 y 2 J d 3 x 3 ^(yi^" 1 (yi, y 2 ) T°' c (x x , x 3 ; y 2 ) 
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G(x 3 ,x 2 ) 



(58) 



r 



where we have again used Wick's theorem Q42p. Finally the defining equation for the ghost-gluon vertex l|51|l we 
expressing the remaining expectation value by means of obtain for the full ghost self-energy 



E(xi,x 2 ) = J d 3 Vl J d 3 y 2 J d 3 x 3 J d 3 x 4 ^f(y 1 ,y 2 )r^(x 1 ,x 3 ;y 2 )G w (x 3 ,x 4 )r£ (x 4 ,x 2 ; yi ) 



(59) 



where we have introduced the short hand notation 



0« (yi.ya) - 2 <fc( ( yi ^ 1 ( yi ' y2 ^ ab ( 6 °) 



for the gluon propagator. Note that and G w are both 
off-diagonal matrices in the adjoint representation. Only 
the gluon propagator (jHDJl, (|4*3|l is color diagonal due 
to our spezific ansatz (1341) . choosing lu color indepen- 
dent. The ghost self-energy S (|59fl is diagrammatically 
illustrated in figure |21 Investigations (in Landau gauge 
|32| ) show that vertex dressing is a subleading effect [36] . 
Therefore we will use in the present paper the so-called 
rainbow-ladder approximation, replacing the full ghost- 
gluon vertex T by its bare one r° (|52fl . Then the ghost 
self-energy becomes 



S ab ( Xl ,x 2 ) = -,5 Qb 7V c3 2 t fei (x 2 V- 1 (x 2 ,x 1 ) 



9?8pG u (xi,x 2 ) 



where we have used 



(f c f c ) 



ah 



= -5 ab N t 



c ■ 



(61) 



(62) 





(a) 



(b) 



FIG. 2: Diagrammatic representation of the ghost self-energy 
eq. 1591 . (a) full, (b) in rainbow-ladder approximation. 



D. The ghost and Coulomb form factor 



Iterating the Dyson equation for the ghost propaga- 
tor we end up with a geometric series in powers of EGo 

00 Id 

Guj = Go } (SGo) = Go- ^77- := Go- . (63) 

^— ' 1 - LG n a 



n=0 



In the last relation we have introduced the ghost form 
factor 



d = 



9 



1-SGo 



(64) 



Using the rainbow ladder approximation l|61|) where £ 
is given by the diagramm shown in figure^b) the Dyson 
equation 149|) for the form factor of the ghost propagator 



becomes 

gd- 1 = 1 - EG := 1 - gh 
or after Fourier transformation 

d(k) g 



(65) 



(66) 



/d(k) - t y (^i 1 -^ J(k-q)Mq)' 

where k = k/fc. Once, the ghost propagator G is known, 
the Coulomb propagator F can be obtained from the 
ghost propagator G by using the relation 



dg 



(67) 



which follows immediately from the definitions J20J), (|17J) . 
This can be rewritten in a more compact form by 



(68) 



Given the structure of the Coulomb propagator it is con- 
venient to introduce yet another form factor, which mea- 
sures the deviation of the Coulomb propagator 



F u := <F% = (G(-a 2 )G), 



(69) 



from the factorized form (G) 0J (—d 2 )(G) 0J . In momentum 
space we define this form factor /(k) by 



F w (k) - G w (k)k 2 /(k)G w (k) 

= 4o d ( k )/( k ) d ( k ) - 
y K 



(70) 



where we have used eq. I|63[) . By taking the expectation 
value of ljrJ5|) we have 



F u = (u\-jL(jgG)\u) 



(71) 



Later on ui will be determined by minimizing the energy. 
Then u> becomes (/-dependent. Ignoring this implicit g- 



dependence of lu we may write (w 
from cq. fTTj). (|7Uj) and |J53J| follows' 



<9G I 
9ff 1 



-G w . Then 



/(k) = ^ d -i (k) ^) d -l (k) = -^^(k) • (72) 

Let us emphasize, that this relation, first obtained in 
|3lj |. is only valid, when the implicit (/-dependence of co 
is ignored. Note also, that in the above equations all 
Greens functions, form factors and vertices like G Ul T,,d 
and / are color matrices in the adjoint representation of 
the gauge group. 

In the rainbow ladder approximation defined by eq. I|61l) 
the ghost self-energy S is given by the one-loop diagramm 
shown in figure fflb) . In this approximation the ghost 



= 1 + 



= 1 




(a) 



(b) 



FIG. 3: (a) Diagrammatic representation of the integral equa- 
tion 1)73(1 for the Coulomb form factor, (b) Series of diagrams 
summed up by the integral equation shown in (a). 



form factor l|66|) is a unit matrix in color space and from 
eq. (|72|) we find for the Coulomb form factor 



/(k) = 1 + 7/Ck), 



(73) 



2 ^ rf (k-q) 2 /(k-q) 
(k-q)Mq) 



Here, we have discarded terms involving |p for reasons 
explained above. The integral equation (|73(l is graphi- 
cally illustrated in figure |3Ja). Iteration of this equation 
yields the diagrammatic serie shown in figure|2Ib). From 
this equation it is also seen, that in leading order the 
Coulomb form factor is given by /(k) = 1. 



E. The curvature 



Consider now the quantity 
5\nJ 6 



<L4£(x) 6A%{x) 



™ ' 1741 



Using the definition of the bare ghost-gluon vertex 15211 , 
we have 



Shi J 



= -Tr 



(GT^(x) 



(75) 



With this relation the momentum operator l|29|l becomes 

n^(x)=n^(x) + iTr(Gr°; Q (x)) . (76) 

Since r° (|54|l is independent of the gauge field, we find 
from 



(77) 



Taking a second functional derivative of eq. I|75|l and 
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(a) 



(b) 





This quantity represents that part of the self-energy of 
the gluon, which is generated by the ghost loop, and, in 
the spririt of many-body theory, is referred to as gluon 
polarization. It can be also interpreted as that part of the 
color dielectric suszeptibility of the Yang-Mills vacuum 
which originates from the presence of the ghost, i.e. from 
the curvature of the orbit space. For this reason we will 
refer to this quantity also as the curvature tensor and the 
quantity 



(c) 



(d) 



FIG. 4: Diagrammatic representation of the curvature. The 
line stands for the full ghost propagator, (a) full curvature 
given in eq. I|78^ . (b) results from (a) by replacing the full 
ghost-gluon vertex V by the bare one I^ ', see eq. I8U1 . (c) 
ghost self-energy correction included in (b). (d) vertex cor- 
rection not included in (b). Note, that all vertex corrections 
contain at least two loops. 



subsequently the expectation value, and using eq. Q51|l. 
we obtain 



5 2 In J 



^(x)Mf(y) 



Tr 

= -Tr 
= -Tr 



ry (y) 



SG 



G w r£(x)G w r°'*(y) 



=: -2x#(x,y) . 



(78) 



x(*,y) 



l l 



2JVg 



1 



5 a6 t«(x)x#(x,y) 



(79) 



is referred to as the scalar curvature |37| . Let us empha- 
size, that the curvature x( x j x ') 1S entirely determined 
by the ghost propagator G(x, x'). Note also, that this 
quantity contains both the full (dressed) and the bare 
ghost-gluon vertices, as illustrated in figure Ufa). 
As already mentioned above from the studies of the 
Schwinger-Dyson equations in Landau gauge [32j | it is 
known, that the effect of the vertex dressing is sublead- 
ing [33. We will therefore replace the full vertex T by 
the bare one l|52|) (rainbow-ladder approximation). The 
curvature is then given by the diagram shown in figure 
S[b). Using the explicit form of the bare ghost-gluon 
vertex (1521 , the curvature tensor (ghost loop part of the 
gluon polarization) becomes 



Xfei(x,y)--ff 2 t fem (x)^(y)Tr 



T a (d*G w (x,y))T fc (c^(y,x)) 



(80) 



r 



In the rainbow ladder approximation used in this paper 
the ghost Greens function is a unit matrix in color space. 
Using tr(f a f b ) = -5 ab N c and x(x,y) = x(x - y) we 
obtain 



t fe „(x) X £?(x,y) = <^i H (x)x(x,y) 



(81) 



with the scalar curvature (|79|l in momentum space given 
by 



d(k-q)rf(q) 
(k - q)2 



(82) 



For the evaluation of the vacuum energy to be carried out 
in the next section we also need the quantity (A 1 - s x ^ ) . 
This expectation value can be conveniently calculated us- 
ing the Wick's theorem in the form of eq. (|42|l and the 



relation (|56l) . We then obtain 
^ a (x')^|y) =\\ ^(^(x'.xi) 



In J) .(83) 



With the definition of the curvature \ (|78|l this quantity 
can be written as 

±a f ^ Sin J 



^r(x') 



<L4£(x) 



d 3 x 1 t fem (x')a;- 1 (x / ) x 1 )x^(xi ) x) . (84) 



We have now all ingredients available to evaluate the vac- 
uum energy. 
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IV. THE VACUUM ENERGY DENSITY 



term) of the Hamilton operator, 
one finds 



Using Wick's theorem, 



For the Gaussian wave functional it is straightforward to 
evaluate the expectation value of the potential (magnetic 



J 



E p = (u\H b \lo) = 



8(0) / d A k 



N C {N ( 



w(k) 



16 



r 



-g 2 5(0) 



d 3 kd 3 k' 



1 



(2?r) 3 u(k)u(k'] 



(3-(kk') 2 ). (85) 



In the above expression for the vacuum energy the di- 
vergent constant 6(0) can be interpreted as the volume 
of space. Removing this constant we obtain the energy 
density, which is the more relevant quantity in quantum 
field theory. 



The evaluation of the expectation value of the kinetic 
term and of the Coulomb term is more involved. Con- 
sider the action of the momentum operator II l)29|l on the 
wave functional 



J 



1 Sin J 



*[A- L ]=iQf(x)*[A- L ], 



where 



<L4£(x) 2<L4£(x) 

Q?(x) = J d 3 xMx,xi)At a (xx) + = J d 3 xMx,xi)At a (*i) - \t ik (*)Tr (gT°'°(x) 



(86) 



(87) 



The kinetic energy (|29|l can then be expressed as 



Analogously the Coulomb energy becomes 



E, 



J d 3 x J d 3 i'(QJ(x)i^«(x)r b (x,x')if' i (x')Q 3 d (x'; 



(89) 



Due to the presence of the curvature term (i.e. the sec- 
ond term) in l|87JI the vacuum expectation value (. 
cannot be worked out in closed form since J [A] is an infi- 
nite series in the gauge potential. For practical purposes 
we have to resort to approximations. Throughout the pa- 
per we shall evaluate propagators to one-loop level and 
accordingly the vacuum energy to two-loop level. To this 
order the following factorization holds in the Coulomb 
term 

(QA ± FA ± Q) U = {F) u [(A ± A ± ) U (QQ) UJ + (AQ) bJ (AQ) UJ ] 

(90) 



where we have used 



(A) u = , (Q) u = 



(91) 



r 



Furthermore in (QQ) U only a single A from one Q has to 
be contracted with an A from the other Q. The remain- 
ing j4's of a Q have to be contracted among themselves. 
All other contractions give rise to diagrams with more 
than two loops. This ensures, that (QQ) U remains a 
total square even after the vacuum expectaion value is 
taken. This fact can be utilized to simplify the evalua- 
tion of (QQ)u 39]. To obtain (QQ) U we need explicitly 
to evaluate only the first term ~ (AA) U and the mixed 
term - 2{A^f£-) u} . The remaining term (^^^f)^ 
can be found by quadratic completion. Using eq. 
we find for the expectation value of the mixed term 
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(^ a (x)Q>(x')> w = ^ a %{x) 



5(x-x')- J d 3 x x uj 1 (x,xi)x(xi,x') 



r 



(92) 



where x(x, x') is the scalar curvature of the gauge orbit by quadratic completion 
space defined by eq. I|79|) . With this result one then finds 



(Q 4 Q (x)Q<;(x')) w = V J d 3 Xl J d 3 x 2 [ W (x,x 1 )-x(x,x 1 )]^(x 1 V- 1 (x 1 ,x 2 )[ W (x 2 ,x')-x(x 2 ,x')] 

The kinetic energy (|88|l becomes then with 5 aa — AT 2 — 1 and t ii (x) = 2 
fiV 2 - 1) 



■ J d3x J d 3x i J d 3x 2 Mx,xi) - x( x 7 x i)] w 1 (xi,x 2 ) [w(x 2 ,x) - x( x 2,x) 



Analogously, one finds for the Coulomb energy using (0, (EJ, (O,® and f abc f abd = N c S cd 
N C (N? - 1) 



E, 



8 



-p 2 / d 3 a; J dVF w (x,x') 



[tyCxJw-^X.x')] 



d 3 xi / d 3 a; 2 [w(x,xi) - x(x,xi)] [tj^x^uj 1 (xi,x 2 )] [w(x 2 , x') - x(x 2 , x')] 

S(x-x')- J d 3 x 1 uj- 1 {x,x 1 ) x (x u x') J Uji(x') 5(x' - x) - J d 3 x 2 u)- 1 (x', x 2 )x(x 2 , x 



iij(x) 



These expressions can be written in a more compact form in momentum space: 



and 



£r = 



TV 2 -, — 1 



<J< 3 >(0) / d 3 k 



'(k)-x(k)f 
w(k) 



WW1,(3) (0) / ^ (l + (kk') 2 ) d(k -^ 2/ J k 2 - k0 



8 



fk-k' 



(Mk)-x(k)]^[^(kQ~x(kQ])' 
w(k)w(k') 



r 



(93) 



(94) 



.(95) 



(96) 



(97) 



where again the divergent factor 5(0) has to be removed functional is determined from the variational principle 

to obtain the corresponding energy densities. 5E(uj)/8lo = 0. Variation of the magnetic energy is 

The kernel w(k) in the Gaussian ansatz of the wave straightforward and yields 

I 



SE P 
5uj(k) 



TV 2 



-5^(0)- 



>(k)2 



k 2 + ^ 

4 y 



d 3 



(#( 3 -^ 2 y w(q) 



(98) 



r 



The kinetic and Coulomb part of the vacuum energy but also implicitly via the ghost form factor <i(k) and 
E c do not only explicitly depend on the kernel w(k), 
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the Coulomb form factor /(k) as well as via the curva- 
ture x(k). However, one can show, that variation of these 
quantities d(k), /(k), x(k) with respect to w(k) gives rise 



to two-loop terms, which are beyond this cope of this pa- 
per. Ignoring these terms, we obtain for the variation of 
the kinetic and Coulomb parts of the energy 



SE k 
<W(k) 
SE C 



%ii^)(o) 



1 - 



X(k) 
W (k)a 



2 1 



8 1 j u(k) 2 7 (2tt) 3 v +l qj 



d(k - q) 2 /(k - q) ^(k) 2 - Mq) - X (q) + X (k)]' 



(k-q) 



;(k) 



(99) 



(100) 



r 



Putting all these ingredients together, we find from the where we have introduced the abbreviations 
variational principle the gap equation to one-loop order 



W (k) 2 = k 2 + X (k) 2 + / W (k) + ^ , 



(101) 



Nc 2 
4 y 



Nc 
4 



d 3 q 
d 3 q 



(2tt) 3 



(kq) : 



1 + (kq) : 



w(q) 

d(k-q) 2 /(k-q) [ W (q) 



(k-q) 2 
I 



X (q)+ X (k)] 2 -c(k) 2 



w(q) 



(102) 
(103) 



We are then left with a set of four coupled Schwinger- 
Dyson equations for the ghost form factor d(k) l|66|l . for 
the curvature x(^) JH2t ; f° r the frequency w(/c) (gap 
equation) l|101|) and for the Coulomb form factor /(fc) 
(|73|l . Before finding the self-consistent solutions to these 
coupled Schwinger-Dyson equations in the next section 
we shall study their analytic properties in the ultraviolet 
(k — > oo) and in the infrared {k — > 0). 



V. ASYMPTOTIC BEHAVIOUR 



of this approximation. However, the numerical calcula- 
tion to be presented in section fVlII which will be carried 
out without resorting to the angular approximation, will 
produce asymptotic behaviours, which are quite similar 
to the one obtained below within the angular approxima- 
tion. 

The angular approximation is defined by approximating a 
function h(k — q), which occurs under the momentum in- 
tegral with an argument given by the difference between 
the external momentum k and the internal loop momen- 
tum q integrated over, by the following expression 



To obtain first insights into the infrared and ultraviolet 
behaviour of solutions of the coupled Schwinger-Dyson 
equation, we investigate these equations in the so-called 
angular approximation, which is defined in eq. (|104f) be- 
low. At this stage we have no estimate of the accuracy 



J 



h(\k - q|) = h(k)G{k -q) + h(q)Q{q - k) . (104) 

With this approximations the angle integrals in the 
Schwinger-Dyson equations for x(^) an d d(k) reduce with 

1 — (kq) 2 ) = sin 2 $ to the following integrals 



7T 

/ 



di9sin d i9 



3 rf(k-q) 



(k - q) 2 



e(k- q m + e(q-k) d -^ 



e(k- q m + e iq -k) d M 



dtfsin 3 i? 



(105) 
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Using these results the remaining integrals over the mod- 
ulus of the momentum occuring in the Schwinger-Dyson 
equations become 



h{k) 



I x (k) = 



N c 
6tt 2 

N c 
12tt 2 



uj(q) 



dq 



d(q) 
uj(q) 



d(k) 



dqq d(q) + / dqd(q) 



(106) 



.(107) 



It turns out, that the simplest way to solve the 
Schwinger-Dyson equation for the ghost propagator and 
the curvature in the asymptotic regions is to differentiate 
these equations with respect to k. The clue is, that the 
derivatives contain only ultraviolet convergent integrals. 
Indeed from equations l|10(j \ , (|107fl we have (The con- 
tributions from the derivative of the integration limits k 
cancel) 



I'M 



W = 



Nc_ _j_ 

6n 2 k 2 

N c 1 
12tt 2 k 2 



d'(k) - 2 
d'{k) - 2 



d(k) 



k 
d(k) 



dq 



Lu{q) 



dqq 2 d{q) .(109) 



Differentiating the Schwinger-Dyson equation for the 
ghost form factor, eq. (|66|l . with respect to the exter- 
nal momentum k and using equation l|108|) . we find 



d'{k) 



1 



N c R(k) 



N c R{k) d(k) 
37r 2 k 2 k 



_d{k) 2 6tt 2 k 2 
where we have introduced the abbreviation 

k 



(110) 



R(k) = J dq 



T 



oj(q) ' 



(111) 



Differentiating the equation for the curvature 1)82(1 and 
using the angular approximation (|107f) we obtain 



X'(k) = 



N c 1 
12tt 2 k 2 



d'(k) - 2 



d(k) 



S(k) 



(112) 



where 



S(k) = j dqq 2 d{q) 




(113) 



Below we will solve these equations separately for k — > oo 
and k — > 0, respectively. 



Ultraviolet behaviour 



Due to asymptotic freedom the gluon energy iv(k) has to 
approach for k — > oo the asymptotic form 



,{k) 



k 2 , k 



(114) 



This behaviour will also be obtained later on from the so- 
lution to the gap equation. Let us therefore concentrate 
on the asymptotic behaviour of the ghost form factor 
d(fc), the curvature x{k) an d t ne Coulomb form factor 
f(k) resorting to the angular approximaton eq. (|105fl . 
Consider first the Schwinger-Dyson-equation (|110fl for 
the ghost form factor. It contains the so far unknown 
integral R(k) 1(111(1 . For large k the dominating contri- 
bution to this integral comes from the large q region. 
Hence, for an estimate of R(k) we can use for uj(k) its 
asymptotic value uj(q) = y 7 ^ 2 yielding 



R(k) = — , k — > oo . 



(115) 



(108) Inserting this result into eq. 1(110(1 we obtain 



d'(k) 



1 



Nr 



d{k) 3 12tt 2 d{k) 



07T Z k 



Integrating this equation yields 



d(k) 



d([i) 



N c A k 2 , d(k)\ , , 

^(v^iy) ' (117) 



where fi is an arbitrary momentum. For large k we ex- 
pect, that the following condition holds 



d(k) k 2 

« — K 



d(fl) 



OO 



(118) 



The reason is, that for k — > oo the Faddeev-Popov deter- 
minant approaches that of the Laplacian and accordingly 
the ghost form factor should approach one, d{k — > oo) — > 
1. However, due to dimensional transmutation giving rise 
to anomalous dimensions, we will shortly see, that d(k) 
behaves asymptotically like d(k) ~ (In k 2 //j, 2 )~ 5 5 which 
of course satisfies the condition eq. 11181) . In anticipation 
of this result, we will in the following assume, that the 
condition eq. 1(118(1 is fullfilled and subsequently show, 
that the resulting solution for d(k) will indeed obey this 
condition. When this condition is fullfilled, we can ig- 
nore the last term in eq. 1117(1 resulting in the explicit 
solution 



d{k) = 



1 



N c 
6tt' 2 



(119) 



which asymptotically behaves like 
d{k) = 



G 



N c 



l &2 



(120) 
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as anticipated before. 

In passing we note, that the obtained asymptotic be- 
haviour of d(k) is precisely that of the running coupling 
constant g ~ ^/a. In section rVfl we will see, that indeed 
asymptotically d(k — > oo) approaches the renormalizcd 
coupling constant. From eq. I|63|) then follows that the 
ghost behaves, indeed, asymptotically like a free particle. 

From the asymptotic form of the ghost form factor H19|) 
we can immediately infer the asymptotic form of the 
Coulomb form factor f(k) by using the relation 
This yields 



/(*) = /(/*) 



d(k) 
W) 



(121) 



Thus, the asymptotic behaviour of the Coulomb form 
factor is up to a numerical constant f (fi) / d(fi) the same 
as the one of the ghost form factor. 

The extraction of the asymptotic form of the curvature 
is somewhat more involved. For large k — > oo we can use 
the asymptotic form (|120[1 of the ghost form factor d{k) 
from which one finds 



d'{k) 



N c 1 
6tt 2 k 



d(k) 3 , 



so that eq. (|112Jl becomes 



X'(k) 



N c d{k) 

67T 2 fc 3 



N c 
12tt 2 



d{kf 



(122) 



S(k) . (123) 



The ghost form factor d(k) is, by definition, strictly pos- 
itive definite inside the first Gribov horizon, to which 
we have to restrict our gauge orbits. Thus the integral 
S(k) (|113JI is positive definite and x'(k) < 0. For suf- 
ficiently large fc, where we can use the asymptotic form 
((HUjl of d{k) for which the integrand in Sik) (|TT3l be- 
haves like q 2 / \Jh\ q/jl. Furthermore, we will show in the 
next section that q 2 d(q) — ► for q — » 40]. Therefore 
for sufficiently large fc an upper limit to S(k) is given by 



S(k) < J dqq 2 = Ifc 3 



(124) 



we obtain the asymptotic estimate 

Integrating this equation over the interval (fco , fc) we find 

> X(k) - x(fco) > -0 N(fc) - k d(k )} . (128) 
Dividing this equation by fc we find that asymptotically 



> 



X{k) 



> 



Nc 
18tt 



or with eq. (tT2U|) 

X(k) 



From this relation it follows, that 

X(k) _ X(k) , 1 
w(fc) fc y/lnk/n 



;d(k) , fc — > oo 



, fc — > oo 



, fc — > oo . 



(129) 



(130) 



(131) 



To summarize, we have found the following ultraviolet 
behaviour (fc — > oo) 



«(k) 

d(k) 

X(k) 



1 



i 



/(k) 



(132) 



The first equation means, that gluons behave asymptot- 
ically (fc — > oo) like free particles with energy fc, while 
the last equation implies, that the space of gauge con- 
nections becomes asymptotically flat. The ghost form 
factor d(k) deviates from that of a free (massless) point 
like particle by the anomalous dimensions 1 / y/lnk/fi. 
These relations are in accord with asymptotic freedom. 

One easily convinces oneself that the asymptotic be- 
haviour obtained above yields indeed a consistent solu- 
tion to the coupled Schwinger-Dyson equation. 



From eq. H123|) we obtain then the following estimate 



0> X '(fc)>-^(fc) 



N c 
12tt 2 



d{kf 



(125) 



Since d{k) ~ 1/ yjlnkj ' \x for fc — > oo the second term in 
the bracket is irrelevant for sufficiently large fc. By the 
same token we can multiply this term by a factor (—2) 
without changing the asymptotic result. Using (|122f> 



d(k) 



1-2^ 



= -N(fc)] 



(126) 



B. The infrared behaviour 



In the following we study the behaviour of the solutions 
of the coupled Schwinger-Dyson equations for fc — » 
thereby using again the angular approximation 11041) . 
For the quantities under interest in the infrared we make 
the following ansatze 



u{k) 



A 



, d(k) 



B 

F 



, X(k) 



C_ 

fc7 



(133) 
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With these ansatze we solve (the derivatives of) the cou- 
pled Schwinger-Dyson equations 1)1101) . 1)112)) for k — ► 0. 
In the remaining integrals 1)111(1 . ()113)1 . the integration 
variable is restricted to the inter vail < q < k. For 
k — > we can use the asymptotic representations, eq. 
I)133)l . in the integrands and obtain 



R(k) = ~ J dqq a+2 = 



S(k) = B / dqq 2 -P = B 



1 fc"+ 3 

Ia + 3 

fc 3-/3 



3-/3 



(134) 



Inserting these expressions into eqs. 1)108)1 . 1)109)1 we find 



6tt 2 A a + 3 



N C „20 + ^ k -2fi 



12tt 2 3-/3 
From the derivative of the ghost equation 

we then obtain the following relation 

A_ = N c /3 + 2 g+1 
B 2 6tt 2 (3(a + 3) 



(136) 



(137) 



The left-hand side of this equation is constant. The same 
has to be true for the right-hand side, which implies 



a = 2/3-1. 



(138) 



Inserting this relation into eq. 1)137)) we obtain a relation 
between the coefficients A and B 



A 



N c (3 + 2 



B 2 6tt 2 2/3(13+ 1) 



(139) 



Analogously we find with 1)135)) from the derivative of the 
curvature equation 



X'(k)=l' a (k) 



the relation 



C N c (3 + 2 



B 2 12tt 2 7(3 - (3) 



(140) 



(141) 



Since the left hand side is a constant (i.e. independent 
of k) it follows 



and thus 



7 = 2/3-1 



C N c (3 + 2 



B 2 12tt 2 (2/3- 1) (3-/3) 



(142) 



(143) 



In view of eqs. 1)138)) and 1)142)1 we have the following 
relations between the infrared exponents 



a = 7 = 2/3 - 1 , 



(144) 



so that ui(k) and x(^) behave in exactly the same way 
in the infrared. Combining eq. 1)139)1 and 1)143)1 we can 
eliminate the constant B and obtain 



A _ (2/3-l)(3-/3) 
C (3((3 + l) 



(145) 



Unfortunately the gap equation 1)101)1 cannot be treated 
in the same fashion. This is because the integrand in 
Iui(k) 1)103)) contains explicitly the external momentum k. 
However, one can show, without resorting to the angular 
approximation, that in the infrared k — > the quantities 
u>(k) and x(fc) approach each other, i.e. u)(k) and x(fc) 
do not only have the same divergent infrared behaviour 
(a = 7) as shown above, but also have the same infrared 
strength 



A = C 



(146) 



For this purpose consider the full gap equation 1)101)) in 
the limit k — > assuming for the moment, that the ultra- 
violet diverging integrals have been regularized. Rcnor- 
malizaton carried out in the next section will remove 
the divergent constant 7° 1)11)2) 1 and will introduce fi- 
nite renormalization constants which, however, become 
irrelevant in the infrared limit k — > compared to the 
diverging quantitites ui(k) and x(fc). Hence, ignoring in- 
frared finite terms, the gap equation ()101)1 becomes 



uj 2 {k -> 0) = X 2 (k -> 0) + /„(* -> 0) 



(147) 



Since the integral 1^ (k) ()103)) is ultraviolet divergent the 
dominant contributions to this integral must come from 
the large momentum region. For large but finite q and 
k — > in the integrant of I u (k) 1)103)1 we can omit x(<?) 
and u>(q) compared to the divergent quantities uj(k — ► 0), 
x(k — > 0) yielding 



J w (fc->0) 



[X 2 (k)~ 

Nc 
4 



A d 2 (k-q)/(k-q) 



27T) 



(k-q) 



(148) 



is ultraviolet finite. Indeed with the above found ultra- 
violet behaviour 1)132)) we obtain 



dq 



d 2 (q)f(q) 
Lu{q) 



dq 



1 



qQnq/fji)* 



dq 



1 



dq (lnq/^)i 
(149) 



With the singular behaviour of u(k) and x(k) for k 
eqs. 1)147)) . 1)148)) obviously imply 



X(k^0) =w(k->Q) 



(150) 
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The same relation is also found in the full numerical so- 
lution of the renormalized gap equation given in section 

IvTTl 

With the relation HZfy it follows from eq. (|TI51) (3 = 1 
and in view of eq. 1144fl a = 7 = 1. Thus, we obtain the 
following infrared behaviour 



uj{k) = X (k) = j , d(k) = 




(151) 



Finally, we note, that the above obtained infrared be- 
haviour for the ghost propagator d(k) ~ i is precisely 
the one, which is needed to produce a linear rising con- 
finement potential from the Coulomb energy, when one 
uses the leading infrared approximation f(k — ► 0) = 1 
for the Coulomb form factor. 



VI. RENORMALIZATION 



The integrals occuring in the Schwinger-Dyson equations 
are divergent and require regularization and renormaliza- 
tion. The renormalization of the Schwinger-Dyson equa- 
tions in Coulomb gauge has been already discussed in 
ref. [23. However, our equations differ from those of ref. 
[23| due to the presence of the curvature x- The latter 
will introduce new features, which require separate dis- 
cussions. In ref. |24j the curvature was also included but 
not fully: In the gap equation the curvature was omitted 
under the loop integral in the Coulomb term, but it is 
precisely this dependence on the curvature, which gives 
rise to the new troublesome features. 

For simplicity, we will use a 3-momentum cutoff A as 
ultraviolet regulator. We are aware of the fact, that 
such a procedure violates gauge invariance and may give 
rise to spurious divergencies. Furthermore, the approxi- 
mate evaluation of the expectation value of the Coulomb 
term (neglecting two-loop terms) will also introduce spu- 
rious ultraviolet divergent terms (see below). The crucial 
point, however, is that neither the infrared nor the ultra- 
violet behaviour of the quantities under interest (deter- 
mined by the above coupled Schwinger-Dyson equations) 
will depend on the specific regularization and renormal- 
ization procedure used as will be demonstrated later. 

After regularization the integrals Id(k), I x (k), I u (k) in 
the Schwinger-Dyson equations become cutoff dependent 
Id(k) — > Id(k,A) etc. . Consider first the equation l|tjtj|) 
for the ghost form factor 



1 



d(k,A) 



1 

fJ 



(152) 



where we have indicated, that after regularization the 
ghost form factor also becomes cutoff dependent. How- 
ever, to assign a physical meaning to d(k, A) we have to 



keep d(k,A) independent of the cutoff A ^l| . This is 
achieved in the standard fashion by letting the coupling 
constant g run with the cutoff A. Independence of the 
ghost form factor d(k, A) of A requires in view of (|152fl 



dg(A) 
dA 



-3 2 (A) 



dI d {k,A) 
dA 



(153) 



For A — > 00 we can ignore the A dependence of u>(q) and 
d(q) in the integrand of Id and find 



dI d (k,A) 



dA 



N c 1 d(fc = A,A) 
3 27r 2 w(fc = A,A) 



(154) 



which is independent of k. Since this quantity is positive 
(oj(k) has to be strictly positive definite and d(k) is pos- 
itive inside the first Gribov horizon) the solution g(A) to 
eq. (I153fl vanishes asymptotically for A — * 00. Further- 
more, for A — > 00 the integral Id{k — A, A) remains 
finite (This is also explicitly seen from the angular ap- 
proximation HID Bp .) Therefore from eq. (|152|) follows, 
that d(k — A, A) approaches asymptotically g(A). Since 
furthermore uj{k — > 00) ~ k we obtain 



A 



dg_ 

'dA 



-Pg 3 , (3 



(4tt) 2 



8N C 



(155) 



This result was also obtained in refs. |23J, [24j, which 
is not surprising, since the curvature does not enter the 
equation for the ghost form factor. As discussed in ref. 
[23 this coefficient should not be compared with the 
canonical perturbative expression of (3q = llNc/3. In 
the present approach the running coupling constant can 
be extracted from the Coulomb term. One finds then 
(ref. 0) /?o = 12N c /3 instead of j3 = UN c /3. The 
difference is due to the absence of the perturbative contri- 
bution due to the emission and absorption of transverse 
gluons, when taking the expectation value of the Hamil- 
tonian. 



Eq. (|155|l has the well-known solution 



5 2 (A) 



/? 5 2( M )m(j£) ' 



1 



(156) 



which shows that for A — > 00 the asymptotic behavour 
of g(A) is 



ff 2 (A) 



(157) 



in accordance with asymptotic freedom. The same be- 
haviour was obtained in eq. I|12l)|) for the ghost form 
factor d(k = A — > 00). This shows, that, indeed, asymp- 
totically the ghost form factor approaches the running 
coupling constant 



d(k = A — > 00) — > g(A) 



(158) 
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The integral Id(k,A) (I66|) is logarithmically ultraviolet 
divergent (This is also explicitly seen in the angular ap- 
proximation p()fi|l .l The equation (|152J) can be renor- 
malized by subtracting the same equation at an arbitrary 
renormalization scale fj, 



1 



d(/i,A) g(A) 



yielding 



1 



d(k,A) din, A) 



[I d (k,A)-I d (fi,A)] 



(159) 



(160) 



Eq. (|159|l shows how (for a fixed renormalization scale fi) 
the coupling constant g(fi, A) has to run with A in order, 
that the ghost form factor d(k) becomes independent of 
the cutoff. For A — * oo this dependence g(A) is given by 
eq. (jT57| . 



is ultraviolet finite, so that the limit A — > oo can be taken. 
The finite quantities x(p) an d f{lA are new renormaliza- 
tion constants, on which our solutions, in principle, will 
depend. However, we will find later, that this dependence 
is very mild and that neither the infrared nor the ultravi- 
olet behaviour of our solutions will depend on x(/x) and 

The renormalization of the gap equation (|l(Jlfl is more 
involved, when the curvature term is included. To renor- 
malize the gap equation we first follow the minimal sub- 
traction procedure (applied above to the ghost form fac- 
tor and to the curvature). This removes the (quadrati- 
cally) divergent constant J° (1 102(1 resulting in 



u{kf 



+Uk) 



k 2 -fi 2 



X(k) 2 - X{»? 



(166) 



The difference Id(k,A) — Id(p, A) is ultraviolet finite, so 
that we can take the limit A — > oo in eq. I|160J1 



1 



1 



d(k) d{n) 



lim [I d (k,A) 

A— >oo 



I d (n, A)] , (161) 



where we have put d(k, A — > oo) = d{k) etc. This is 
the desired finite Schwinger-Dyson equation for the ghost 
form factor, which contains the arbitrary renomalization 
constant d{[i). 

We use the same minimal subtraction procedure to renor- 
malize the equation for the curvature (|82|) and the 
Coulomb form factor i|73|) . Subtracting the equation 182|) 
and H73f) once at the renormalization scale [i yields 



x(*)=x(/0+A/ x (fc) 
f(k) = f( f i) + AI f (k) ) 



where the difference 



(162) 
(163) 



Unfortunately, the resulting expression 7 w (fc) — Iu{[J.) is 
still diverging. This is a consequence of the one-loop 
approximation, used when calculating the expectation 
value of the Coulomb term (resulting in the factoriza- 
tion (j^UJl) and also when taking the variation (functional 
derivative) of this term with respect to ui(k) to obtain 
the gap equation. As is well known the gauge invariance 
is not maintained in the loop expansion order by order 
|27| . As a consequence, truncating the loop expansion at 
a given order results in spurious divergencies, which are 
cancelled by higher order terms. Such spurious divergent 
terms should therefore be omitted. To identify the spu- 
rious (divergent) terms in the gap equation, we rewrite 
the diverging integral in the form 



Uk,A)=lW{k,A)+2x(k)I^{k,A) 



(167) 



A7 x (fc) = 7 x (fe,A)-7 x (/i,A) 
AI f (k)=I f (k,A)-I f (»,A) 



(164) 
(165) 



where 



lL n) (k,A) = 



N c 



d 3 q 
(2nY- 



(i + (kq) 2 ) 



d(k - q) 2 /(k - q) Mq) - x(q)]" - [w(k) - X (k)f 



(k-q) 2 



w(q) 



(168) 



r 



(2) 

The integral 1^ (fc) is quadratically divergent. Its diver- 
gent part is, however, independent of the external mo- 
mentum fc, so that one subtraction eliminates this diver- 
gence, i.e. 1^ (k) — 1^ ' (/i) is finite. The troublesome 



term in eq. (|167f) is the second one, which is linearly 
divergent. Due to the momentum dependent factor x(^) 
one substraction does not eliminate the divergency 
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J w (fc,A)-/ w ( M ,A) 



I™ (fc, A) - /( 2 ) (,;, A)] +2 x (fc) [l™ (fc, A) - ( M , A)] +2 [ x (fc) - X (/x)] ( M , A) , (169) 

divergent 



finite 



finite 



l!u\ix) is finite, the last term is still diverg- which we denote by l't} x> (fx). In the following we will 

only keep this finite part i!} (fx), which unfortunately is 
not explicitly known. However, we will later show, that 
this unknown constant l'S^ (/i) does not essentially influ- 
ence the solutions of the renormalized coupled Schwinger- 
Dyson equations. 

Applying the above described renormalization procedure 
(i.e. one subtraction and replacing l£P (/i, A) by its finite 
part l'u r (m)) the g a P equation l|l()6|l becomes 



while I { J } (k)- 

ing. Note, that this term disappears, when the curvature 
is ignored (x(fc) — > 0) as done in ref. 23] 42]. (How- 
ever, we will later find that it is of crucial importance 
to fully keep the curvature.) Due to its fc dependence a 
further substraction would not eliminate this divergency. 
This is the type of spurious term discussed above, whose 
singular part violates gauge invariance and would be can- 
celled by higher order loop terms. Once these higher or- 
der loop terms are included the divergencies are cancelled 
and what is left from I^\/x,A) is a finite contribution, 



uo(kf = k 2 + X {kf+uo( j xf-^- X ^f + [/(, 2 >(fc,A)-4 2) taA) 
+2 x (fc) [/« (fc, A) - /« (//, A)] + 2 [ X (k) - X (/z)] I'V> (fx) . 

Inserting here for X (fc) its renormalized value (|162f) we obtain 

c(fc) 2 = fc 2 - A1 2 + A/ x (fc) 2 +eA/ x (fc)+[/( 2 )(fc,A)-/( 2 )( M ,A) 
+2[ X (n)+AI x (k)} f/W(fc, A) -!«(/,, A)l + ^) 2 . 



(170) 



(171) 



All together there are five renormalization constants 

<2(m),x(/"0,/G"),w(/z) and £ = 2[x(/x) + (/*)]• Later 
one we shall demonstrate that the self-consistent solu- 
tion to the coupled Schwinger-Dyson equations does not 
sensitively depend on the detailed values of these renor- 
malization constants except for d(fi). 
From our previous discussions it should be clear, that the 
ultraviolet behaviour of the self-consistent solution does 
not at all depend on these renormalization constants. We 
will later also show, that the infrared behaviour does not 
depend on the precise value of the renormalization con- 
stants x(a*)j /(/•*), ^(a 4 ) an d while the ghost form factor 
d(k) depends crucially on d(fx) and only for one particular 
(critical) value of d(/x) (for which l/d(fx — > 0) — * 0) the 
coupled Schwinger-Dyson equations will have a solution, 
consistent with the infrared behaviour found in section 
IVBl 



VII. NUMERICAL RESULTS 



In this section we present the results of the numerical so- 
lutions to the coupled Schwinger-Dyson equations (|161fl . 
dTTTJ, (UH2), fH%|l for the ghost form factor d(k), the 
gluon energy w(fc), the curvature X (k) an d the Coulomb 
form factor /(fc). For this purpose it is convenient to 
introduce dimensionless quantities. We will rescale all 
dimensionfull quantities with appropriate powers of the 
gluon energy 6 :— u>(fx) at an arbitrary renormalization 
point fx. The rescaled dimensionless quantities will be 
indicated by a bar 

- fc _ u(k = kS) - x (k = kS) 

k = g , u(k) = , XW = g ■ (172) 

The form factors d(k) and /(fc) are dimensionless. 
Before solving the coupled Schwinger-Dyson equations 
we have to fix the renormalization constants 

d(fx) , £ = 2[ x (fx) + I'M (fx)} , x(M) , m • (173) 

Note, that ui([x) has been absorbed into the dimension- 
less quantities. 
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FIG. 5: Solution for the ghost form function d(k) for different 
renormalization constants d(n) = 7.0, 8.6 and 8.716. 



The coupled Schwinger-Dyson equations are solved by 
iteration without resorting to the angular approxima- 
tion. To carry out the calculation consistently to one 
loop order, we should use the leading expression for the 
Coulomb form factor f(k) = 1 (see eq. fZHJ)- How- 
ever, then we loose the anomalous dimension of /(fc)(~ 
l/-\/ln k 2 //j, 2 , k — > oo) which is needed for the conver- 
gence of certain loop integrals. To keep the anomalous 
dimensions of f(k) and at the same time include as lit- 
tle as possible corrections to the leading order f(k) = 1 
we replace the ghost form factor d(k) in the equation for 
f(k) by its bare value d(k) = 1 and use /(//) = 1 . 
The integrals were calculated by using the Gauss- 
Legendre method. In order to obtain an accurate map- 
ping of the infrared region, a logarithmical distribution 
of the supporting points was used. The self-consistent 
solutions are shown in figures EHHI and [7| for the choice of 
the renormalization constants £ = 0, x(/i) = 0. 
The value of the remaining renormalization constant 
d{n) has been specified as follows: 

Consider the equation for the ghost form factor. The 
curvature x(fc), the Coulomb form factor f(k) as well as 
the renormalization constants £ and x(m) do not enter 
this equation. Thus, for given Q(k) the solution d(k) de- 
pends only on the renormalization constant d(/i). Figure 
[5] shows the solution to the ghost form factor for various 
values of the renormalization constant d(fi) keeping ui(k) 
fixed to the solution shown in figure |SJ It is seen, that all 
solutions have the same ultraviolet behaviour indepen- 
dent of the renormalization constant d{k). Furthermore, 
this ultraviolet behaviour is consistent with the asymp- 
totic solutions l|120f> found in section (Note, that a 
double logarithmic plot is used) ! The infrared behaviour 
of d(k) depends, however, on the actual value of d(fi). 
For d(n) smaller than some critical value d cr the curves 
approach a constant for k — > 0. At a critical d(fi) = d cr 
the ghost form factor diverges for k — > and above the 
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Solution for the gap function u>(k) for £ = 




FIG. 7: Ghost form function d(k) and Coulomb correction 
f(k) for£ = 



critical value d > d cr no solution to the ghost form factor 
exists. We have adopted the critical value d(fi) — d cr as 
the physical value for the following reasons: 

i) In D = 3 (which will be considered elsewhere) a 
self-consistent solution to the coupled Schwinger- 
Dyson equations exists only for this critical value. 

ii) Only the critical value produces an infrared diverg- 
ing ghost form factor. 

iii) The diverging ghost form factor is in agreement 
with our analytic studies of the infrared limit of the 
Schwinger-Dyson equations considered in sectionlVl 
using the angular approximation (see eqs. I|151|) ). 

iv) The divergent ghost form factor gives rise to a lin- 
ear rising confining potential as will be shown later 
on. 

The critical d(fi) = d cr is defined by d _1 (fc — > 0) — > 
which is referred to as "horizon condition" |2a|. At the 
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FIG. 8: Ghost form function d(k) for f = 0, -0.5 and -1.0 . 
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FIG. 9: Gap function Q(k) for f = 0, -0.5 and -1.0 . 



arbitrarily chosen (dimensionless) renormalization point 
p, = 0.32 the critical renormalization constant is given by 
d(fi) = 8.716. In all self-consistent solutions presented in 
this paper we have adopted this critical value. 
We have also investigated the dependence of the self- 
consistent solutions on the remaining renormalization 
constants £ and x(/Lt) (recall, that d,u>(k) are indepen- 
dent of £ and xin))- We have found, that our self- 
consistent solutions change by less then 0.01°/oo, when 
x(/i) is varied in the intcrvall [—1, 1]. Thus there is prac- 
tically no dependence of our results on x(fJ,)- We have 
therefore put x(/z) = in all calculations. 
Figures |H1 and [5] show the self-consistent solution for d(k) 
and Q(k) for £ = 0,-0.5,-1.0. Both quantities show 
only very slight variations with £ up to a (dimension- 
less) momentum of order one. The ultraviolet behaviour 
is independent of £ and in agreement with our analytic 
results obtained in section [V] Furthermore, also the in- 
frared behaviour of d(k) and ui(k) is independent of £. 

Our analysis of the infrared behaviour of the solutions 
to the Schwinger-Dyson equations using the angular ap- 
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FIG. 10: A/B 2 and a in dependence of f 



proximation given in section Ivl has revealed, that in this 
approximation the critical exponents a and (3 and the 
ratio A/B 2 of the amplitudes of ui(k) and d(k), see eq. 
(|133[1 . are independent of the renomalization constants 
£, x(/i). Our numerical solutions confirm this result even 
without resorting to the angular approximation. Figure 
1101 shows these quantities A/B 2 and the infrared expo- 
nent (|133[) a as function of £. There is practically no 
dependence. We will therefore put £ = in the further 
numerical calculations. 

The obtained numerical results are all in qualitative 
agreement with our previous analytic investigations. For 
large k — > oo the gluon energy uj{k) ~ Vk 2 is that of a 
non-interacting boson and the curvature in orbit space 
x(k)/u>(k) ~ —t 1 vanishes asymptotically. This is in 

agreement with the expectations of asymptotic freedom. 
For k — > the gluon energy cu(k) diverges reflecting the 
absence of free gluons in the infrared, which is a man- 
ifestation of confinement. While the gluon propagator 
— — > for k — > is suppressed in the infrared the 

ghost propagator d(k)/k 2 diverges for k — > 0. It is worth- 
wile noticing, that the same behaviour of the gluon and 
ghost propagators is obtained in covariant Schwinger- 
Dyson equations, derived from the functional integral in 
Landau gauge |23], |29|- From figure it is seen, that 
for k — > u>(k) approaches x(k). As shown analytically 
in section this is a generic feature of our gap equa- 
tion and reflects the non-trivial metric of the space of 
gauge orbits (given by the Faddeev-Popov matrix). This 
non-trivial metric is crucial for the infrared behaviour of 
the theory and in particular for the confinement. This 
can be seen from figures ^2 1121 > where we present the 
self-consistent solutions for u>(k) and d(k), when the cur- 
vature of the gauge orbit space is neglected by putting 
x(k) = as done in 23] or when neglecting the curvature 
in the Coulomb term as done in |2J]. In these cases the 
infrared behaviour of u>(k) is drastically different from the 
previous case, although we have still chosen the horizon 
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FIG. 11: Consistent solution of the gap function Q(k) for 
different treatments of the curvature for £ = . 
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FIG. 12: Consistent solution of the ghost form function d(k) 
for different treatments of the curvature for £ = . 

condition d~ 1 (k — > 0) = as renormalization condition 
(d(k) is still infrared divergent as can be read off from 
figure IT2*|) . In particular notice that oj(k — * 0) = const, 
when the curvature \ is completely neglected, i.e. a = 
in eq. 1133|) . From the two sum rules 1)1441) for the in- 
frared critical exponents follows (5 = \ and 7 = 0. Thus 
with the horizon condition as renormalization the neglect 
of the curvature in the Schwinger-Dyson equation yields 
d(k) ~ -j=, co(k) = const, and xify = const, for k — > . 



VIII. THE COULOMB POTENTIAL 



The vacuum expectation value of the Coulomb term of 
the Yang-Mills Hamiltonian can be interpreted as in- 
teraction potential between static color charge densities 
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FIG. 13: Coulomb Potential for f = with (full line) and 
without inclusion of the curvature (dashed line). 

p a (x). The static quark potential can therefore be ex- 
tracted from this term by taking the vacuum expectation 
value and assuming, that the color charge density p°(x) 
describes two static infinitely heavy color charges 

P a {x) = gtfrf (x - x (1) ) + gqf 2) 5 (x - x (2) ) (174) 

located at x^) and X( 2 ) and separated a distance xm — 
X( 2 ) = r apart. This yields 

E C - 4 1) +4 2) +^ 1) y Qb (x (1) ,x (2) )^ 2) ,(175) 

where E y c ' ' are the (divergent) self-energies of the two 
separate static quarks and 

V ab (x (1) ,x (2) ) = g 2 (La\F ab (x (1) ,x (2) ) | W ) (176) 

is the static quark potential, with F ab (x,x') being the 
Coulomb propagator defined by eq. fT7| . In the above 
considered one-loop approximation the potential is color 
diagonal V ab = 6 ab V and, with the explicit form of the 
Coulomb propagator, is given by 



V(r) = 




V{k).(177) 



Performing the integral over the polar angle, one finds 

CO 

nr) = ±J dkd{kf m ^ . (178) 



Before presenting the numerical result for the Coulomb 
potential let us consider its asymptotic behaviour for 
k — > and k — > oo. In section IV Bl we have found the 
infrared behaviour f(k — * 0) = const, and d(k — > 0) ~ p 
This yields precisely a linearly rising Coulomb poten- 
tial V(k) r~j 1/fc 4 . Furthermore for k — > oo the ghost 
form factor was found to behave as (see eq. I|12(J|) 
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d(k) ~ 1 Adopting the leading order expression 

for the Coulomb form factor (see eq. (|73J) and figure OH) 
f{k) = 1 we find 

k 2 V{k) ~ — !— , fc -> oo . (179) 
In fc/^ 

This is precisely the behaviour found in ref. [30l | in one- 
loop perturbation theory. 

The Coulomb potential calculated from the numerical 
solution to the coupled Schwinger-Dyson equations is 
shown in figure ^] At small distance it is dominated 
by an ordinary ~ ~ potential, while at large distances 
it raises almost linearly. The numerical analysis shows, 
that its Fourier transform behaves for k — > like 1/fc 3 7 , 
while a strictly linearising potential would require a l//c 4 
dependence (In ref. [34J the power l/k 3 - 6 was found). 
When the curvature is neglected the gluon energy be- 
comes infrared finite and the Coulomb potential ap- 
proaches a constant at r — > oo. Thus both quark and 
gluon confinement is lost when the curvature of the space 
of gauge orbits is discarded. 



IX. SUMMARY AND CONCLUSIONS 



and gives rise to a linear rising static quark potential. 
The asymptotic analytic solutions for both k — ► and 
k — ► oo are reasonably well reproduced by the full nu- 
merical solutions of the coupled Schwinger-Dyson equa- 
tions. Our investigations show, that the inclusion of the 
curvature, i.e. the proper metric of orbit space, given 
by the Faddeev-Popov determinant is crucial in order to 
obtain the confinement properties of the theory. When 
the curvature is discarded (using a flat space of gauge 
connnections) free gluons exist even for k — ► and the 
static quark potential is no longer confining. 

The results obtained in the present paper are quite en- 
couraging and call for further studies. In a subsequent 
paper we will investigate along the same lines the 2 + 1- 
dimensional Yang-Mills theory, which (up to a Higgs 
field) can be considered as the high temperature limit 
of the 3 + 1-dimensional theory. It would be also in- 
teresting to calculate the spatial Wilson loop in order to 
check whether the relation a cou i ~ 3cr found on the lattice 
[lH is reproduced. Furthermore, the spatial t'Hoft loop 
should be calculated using the continuum representation 
derived in |33| . Eventually one should include dynamical 
quarks, since the ultimate goal should be the description 
of the physical hadrons. 



In this paper we have solved the Yang-Mills Schrodinger 
equation for the vacuum in Coulomb gauge by the varia- 
tional principle using a trial wave function for the Yang- 
Mills vacuum, which is strongly peaked at the Gribov 
horizon. Such a wave functional is recommended by the 
fact, that the dominant infrared field configurations lie 
on the Gribov horizon. Such field configurations include, 
in particular, the center vortices, which have been iden- 
tified as the confiner of the theory. With this trial wave 
function the vacuum energy has been evaluated to one- 
loop order. Minimization of the vacuum energy has led 
to a system of coupled Schwinger-Dyson equations for 
the gluon energy, the ghost and Coulomb form factor 
and for the curvature in orbit space. Using the angu- 
lar approximation these Schwinger-Dyson equations have 
been solved analytically in both, the infrared and the ul- 
traviolet regime. In the latter case, we have found the 
familiar perturbative asymptotic behaviours. In the in- 
frared the gluon energy diverges indicating the absence of 
free gluons at low energies, which is a manisfestation of 
confinement. The ghost form factor is infrared diverging 



Note added 



After this work was completed we have been able to show 
that the infrared limit of the Yang-Mills wave functional 
(|3u[) is independent of the power a of the Faddeev-Popov 
determinant |35|. 
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